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We model the secular evolution of a star's orbit when it has a nearby binary 
system. We assume a hierarchical triple system where the inter-binary distance is 
small in comparison with the distance to the star. We show that the major secular 
effect is precession of the star's orbit around the binary system's centre of mass. We 
explain how we can obtain this precession rate from the star's radial velocity data, 
and thus infer the binary system's parameters. We show that the secular effect of 
a nearby binary system on the star's radial velocity can sometimes mimic a planet. 
We analyze the radial velocity data for t/-octantis A which has a nearby companion 
(i/-octantis B) and we obtain retrograde precession of (—0.86 ± 0.02)°/yr. We show 
that if i/-octantis B was itself a double star, it could mimic a signal with similarities 
to that previously identified as a planet of t^-octantis A. Nevertheless, we need more 
observations in order to decide in favor of the double star hypothesis. 



1 INTRODUCTION 



Most extra-solar planet detections rely on measuring the 
parent star's wobble which is assumed to be caused by a 
planet. However, other effects can cause stellar wobble thus 
it is important to study these in order to avoid erroneous 
new planet announcements, as it already h appened in the 
past (jOueloz et al.ll200ll : ISantos et allbOO^ ) 

In previous articles i Morais fc Correiall2008l . l201ll 'l. we 
studied the short-term effect of a binary system on a star's 
motion. We saw that this could mimic a planet companion 
to the star under some circumstances. In these articles, we 
considered moderately close binary systems (> 10 AU) in 
which the star's motion around the binary's centre of mass 
had periods of several decades. Moreover, we realistically 
assumed that we had observational data for only a fraction 
of this period, and not several orbits. 

Here, we will present another scenario that requires a 
different analysis. We consider a star with a close binary 
system (< 5 AU) and we assume that the observational data 
covers a few periods of the star's motion around the binary 
system's centre of mass. We will show that, in this case, we 
have to take into account secular effects which lead to slow 
precession of the star's orbit. 

In Sect. 2 we present the secular theory for hierarchical 
triple star systems composed of a star and a nearby binary. 
In Sect. 3 we show how we can measure the secular preces- 
sion of the star's orbit from radial velocity data and how we 
can predict the binary system's parameters from this mea- 
surement. In Sect. 4 we apply the results from previous Sec- 
tions to fictitious hierarchical triple star systems. In Sect. 5 
we discuss the reported finding of a planet in the binary 
system //-Octantis. In Sect. 6 we present our conclusions. 



2 SECULAR THEORY FOR HIERARCHICAL 
TRIPLE STAR SYSTEMS 

We consider a triple star system composed of an observed 
star, m2, and a nearby binary of masses mo and mi. We 
use the Jacobi coordinates rl (distance of mi to mo), and 
(distance of m2 to the centre of mass of mo and mi). 
Moreover, we assume that |ri| ^ |r2| (hierarchical triple 
system) . 

2.1 Secular Hamiltonian 



The Hamiltonian is (|Lee fc Peald I2OO3I : iFarago fc Laskaj 
l20ld ) 



H = 



Gmorrii G(mo -f- mi)m2 



+ F 



(1) 



2ai 2 02 

where the 1st term describes the Keplerian motion of mi 
with respect to mo (inner binary), the 2nd term describes 
the Keplerian motion of m2 with respect to the centre of 
mass of mo and mi (outer binary), and 

F = — Gmom2 ( — — 

\ ro2 r-2 , 

where G is the gravitational constant, the distance of m2 to 
mo is 

mi 



Gmim2 



Vri2 r2 



(2) 



ro2 =r2 + 



-ri 



mo -I- mi 

and the distance of m2 to mi is 
mo 



ri2 — r2 



-ri 



(3) 



(4) 



mo + mi 

Expanding l/ro2 and l/r-i2 in powers of p = '"i/'"2, and 
retaining terms up to order p^, we obtain the quadrupole 
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Hamiltonian 

G'm.2 momi 



AG 



F 



r2 mo + mi 2 



(3 (fi-fa)' 



(5) 



where fi and f 2 are the versors of fi and r2 , respectively. 

The secular quadrupole Hamiltonian is obtained by av- 
eraging Eq. ([SI with respect to the inn er and outer binary's 
orbital periods (|Farago fc Laskarll2010l ) 



F 



where 



12e?-6(l-e?)(fci-fe2)^ 



c [2 

+30ei{ii-k2f 



C = 



G m^mi 



m2 



16 mo + mi (1 



ei)3/2 al ' 



(6) 



(7) 



fc2 and fci are, respectively, the versors of the angular mo- 
mentum vectors of outer binary (G2) and inner binary (Gi), 
and i\ is the unit vector in the inner binary's orbital plane 
that points towards the inner binary's pericentre. 
In an arbitrary reference frame we have 

fci • fc2 = sin/i sin/2 cos(ni — 0,2) + cos/i COSI2 (8) 
ii ■ k2 = — sinJ2 cos/i sinoji cos(ni — fl2) 

— siriJ2 costal sin(f2i — fl2) 

+ sin Ji sin uji cos I2 , (9) 

with Ii (inner binary's inclination), I2 (outer binary's in- 
clination), Qi (inner binary's longitude of ascending node), 
il2 (outer binary's longitude of ascending node), lji (inner 

bin ary's ar gu ment o f per icentre) . 

Kozail ildei ): iKrvmolowski fc Mazehl (|l999h : 



I Ford et all (|200Cll ^ write the Hamiltonian in the invari 
ant plane reference frame (Fig. 1). In this reference frame, 
^li ~ i^2 — 180°, and the relative inclination is i = Ii + I2, 
thus 

F = C'[(2 3ei)(3 cos^ j - 1) -I- 15ei sin^ i cos(2 wi)] (10) 
and due to conservation of angular momentum 

const . (11) 



= GI+GI + 2Gi G2 cosi 



2.2 Secular equations 

The Delaunay canonical variables for this triple system are 
the angles Ij (mean anomalies), Uj (arguments of pericen- 
tre), Qj (longitudes of ascending nodes), and their conjugate 
momenta, respectively 



Gi = L,J\-e 



with j = 1 (inner binary) and j — 2 (outer binary), 
momi 



mo -I- mi 
m2(mo + mi) 



(12) 
(13) 
(14) 

(15) 
(16) 



Pi = 
P2 = 

mo + mi + m2 

fii — G(mo + mi) and fi2 = G(mo -I- mi -I- m2). 

By definition, the secular Hamiltonian does not depend 
on the mean anomalies, Ij, hence their conjugate momenta. 




Figure 1. The invariant plane is orthogonal to the total angu- 
lar momentum vector G = Gi + G2, where G2 and Gi are, re- 
spectively, the angular momentum vectors of the outer and inner 
binary's. 



Lj, and thus the semi-major axes, aj, are constant. The 



secular evolution is obtained from 



OF 

dujj 
OF 
dG, 

dF 

an^ ' 

OF 
"dlF 



(17) 



(18) 



(19) 



(20) 



The quadrupole Hamiltonian (Eq. [6]) does not depend on ijJ2 
hen ce from Eq. (|17p with ,7 = 2, G2 and 62 are constant. 

iKozail (|l962f l and lKinoshita fc Nakail (|2007l ) derived ex- 
pressions for the secular evolution of the inner binary in 
the limit mi = (inner restricted problem), using the 
quadrupole Hamiltonian. In this case, the outer binary's or- 
bit coincides exactly with the invariant plane which is the 
natural choice of reference frame. The Hamiltonian is given 
by Eq. HlOfFl and the secular evolution of the inner binary i.e. 
{ei,uui) are given by Eqs. (|17p and 1)18^ with j = 1. More- 
over, the secular oscillations of ei and i are coupled due to 
conservation o f angular mo ment um fEq. | ll[p | . 

Following iKozail l|l962h and iKinoshita fc Nakail (|2007l ) 
we will describe the secular dynamics of prograde orbits 
(0° < i < 90°) but since Eq. (|10|l is invariant with respect 
to the transformation i — > 180° — i, the secular dynamics is 
the same for prograde (inclination i) or retrograde (inclina- 
tion 180° — i) orbits. When < i < ic ~ 40°, oJi circulates 



^ When mi = we must replace /3i = 1 (Eq. I15t into Eq. [7] 
llKozailll962l). Note th a t ther e is a mistake in the expression given 
in IKinoshita &: Nakail ||2007| ) since it should have a| and not 
in the nominator. 

^ Since G^ = Gf + G2 + 2 Gi G2 cos i = const, G2 = const and 
Gi ■C G2 hence (1 — ej) cos^ i const 
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while the eccentricity, ei, and relative inclination, i, exhibit 
secular oscillations with amplitude that increases with the 
relative inclination (in particular, coplanar orbits keep con- 
stant eccentricity, ei). When i > ic~ 40°, there are station- 
ary solutions (i — const, ei = -y/l — (5/3) cos^ i — const, 
wi = ±90°), and Kozai cycles where i, e\ and oscillate 
aro und the stationary solutions. 

iFarago fc Laskaij (|2010l ) showed that the secular motion 
of the inner binary when < i < 90° (prograde orbits) is, 
in the general problem (mi 7^ 0), equivalent to the inner 
restricted problem (mi =0), as long as 

G?-L?+4G2 + 2Gi G2 cosi>0 . 



(21) 



Defining X = L2/I/1, we can write Equation (I2ip as 



4(1-62) X'^ + 2yjl - ei \J\-el cosiX -el > . (22) 



The left-hand side of Eq. (|22p is a second degree polyno- 
mial in X which is a convex function with two roots Xi < 
and 5/2 > X2 > Xi. Therefore, inequality (Eq. [511 is ver- 
ified if L1/L2 < 2/5 which is generally true if 01/02 <^ 1 
(hierarchical system) unless mo + mi 3> m2 . 

Equations (|17|l and (|18|l with j = 2 describe the secular 
evolution of the outer binary's eccentricity and argument of 
pericentre. We saw that the secular quadrupole Hamiltonian 
(Eq.O does not depend on (jJ2, hence from Eq. (I17p with j — 
2, the conjugate momentum, G2, and thus the eccentricity, 
62, are constant. From Eq. (llSp with j = 2 we obtain the 
outer binary pericentre's precession rate 



UJ2 



12 c 

15 



1 
2 

2 f~ 

6i(li 



o 2 
d 61 



k2 



(1- 



e?)(fci 



k2 



(23) 



Higher order secular octupole terms cause long-term 
small amplitude oscillations in 61 and 62, which are more 
important f or small to moderate values of the relative in- 
clination, i jKrvmolowski fc Mazehl Il999l : iFord et all I2OO0I : 
iLee fc Pealell2003l '). 



2.3 Precession of the outer binary's orbit 

The outer binary's precession rate (Eq. I23p depends on the 
secular motion of the inner binary. Moreover, in the invariant 
reference frame (Fig. 1) we have 



Gi sin I\ = G2 sin I2 



(24) 



Since a\ <gi 02 (hierarchical system) then, in general, Gi ^ 
G2 (unless mo -I- mi ^ m2) thus sin/2 ^ 1 i.e. the outer 
binary's motion coincides approximately with the invariant 
plane. Therefore, we express the right hand side of Eq. (|23p 
using the reference plane of the outer binary's orbit, i.e. 
setting 72 = and Ji = i in Eqs. ([8}, (|9]) thus obtaining 



UJ2 



A 



12 G 

"g7 



A 



(i.|ef)(3, 



1) 



(25) 
(26) 



+ i^6?(l-e2)cos(2c.i) , 



where 9 = cos i. 

Equation (|25p is an approximation of the precession 
rate, (1)2, because the outer binary's orbit is not fixed but 



exhibits small amplitude oscillations around the invariant 
plane. The angle u)\ on the right hand side of Eq. (|26p 
is measured with respect to the outer binary's orbit or, 
equivalently, with respect to the invariant plan^. This 
formulation is neces sary in order to describe the motion 
of the inner binarv (iKozail Il962l : iKinoshita fc Nakail l2007l : 
iFarago fc Laskaij I2OI0I ). However, the angle 0^2 represents 
the location of the outer binary's periapse with respect to 
the intersection with the observer's plane (when dealing with 
radial velocity data, this is the plane orthogonal to the line 
of sight). 

The long-term evolution of Kozai cyles 
(which exist i f ^ > ~ 40°) w as i nvest i- 

gated bv lEggleton fc Kiseleva-EggletonI (|200lh : 
iFabrvckv fc Tremaind (l2007l ): IW u et all (l2007l ). Typi- 



cally, if 61 becomes close to unity during a Kozai cycle, a 
combination of tidal evolution and relativistic effects will 
eventually disrupt the Kozai cycle and freeze the relative 
inclination, i. This will be followed by tidal damping of the 
semi-major axis, a\, and eccentricity, ei. The end state of 
a Kozai cycle that reaches 61 « 1 will be a tighter inner 
binary (smaller a — 01/02) on a circular orbit. Obviously, 
if Q 1 then, as dj2 oc n2, the precession of the outer 
binary's orbit will be slow thus difficult to detect from 
observational data. On the other hand, orbits nearby the 
Kozai stationary solution are less prone to undergo tidal 
evolution and should keep the original value of a. 

We will, therefore, assume three scenarios for the inner 
binary's motion in the invariant plane reference frame : 

• i < ic where ei and 9 are average values of the secular 
oscillations in ei and 9, respectively. If ei 7^ then uji circu- 
lates and we have, on average, A = (1/2-1-3 ei/4)(3 9^ — 1). 
Uei^O then A « (3 - 1) /2. 

• i > ic but the inner binary's orbit was initially a high 
amplitude Kozai cycle that was circularized by tidal damp- 
ing. In this case we also have A= {3 9^ — l)/2. 

• i > ic and the inner binary is at the Kozai stationary 
solution with loi — ±90° and 9^ — 3{1 — ei)/5. In this case 
A = -5{29^ - 1){9^ - 1). 

In all scenarios above the precession rate, (1)2 is approxi- 
mately constant. 

In Fig. 2 we plot the normalized precession rate, A given 
by Eq. (|26p . when ei = and at the Kozai stationary solu- 
tion. We see that when ei = 0, precession is prograde when 
i < 54.73° and retrograde when i > 54.73°. At the Kozai 
stationary solution, which exists only when i > ic ^ 40°, 
precession is retrograde when i > 45°. 

Figure 2 shows the normalized precession rate. A, for 
0° < i < 90° (prograde orbits). However, we saw previously 
that Eq. (|10p is invariant with respect to the transformation 
i — ^ 180° —i, hence the precession rate (Eq. refatheta) is the 
same for prograde (inclination i) or retrograde (inclination 
180° - i) orbits. 

We performed numerical integrations of the equations 
of motion of hierarchical triple star systems with parameters 



^ Since the intersection of inner and outer binary's orbits (line of 
nodes) is in the invariant plane (Fig. 1), then the angle wi is the 
same when measured with respect to the outer binary's orbit or 
with respect to the invariant plane. 
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■ ■ ■ ■ ' Kozai solution 


□ a=0.1 
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a=0.15 


o a=0.15 



Figure 2. Comparison between theoretical precession rates {A 
given by Eq. I26|l and values obtained in simulations. 



m2 — Mq, a2 — 3.0 AU, 62 = 0.2, and initial angles I2 ~ 0° , 
0,1 = fl2 and cJi = 90°. The inner binary had mi = 0.08 Mq 
and mo = 0.42 M©, semi- major axis ratio a = ai/02. We 
chose two configurations: (i) i = and ei = 0.01 (coplanar 
nearly circular orbit); (ii) i — 60° and ei — 0.76 (Kozai sta- 
tionary solution). In Fig. 2 we show the comparison between 
the theoretical precession rates and the values obtained in 
the simulations (up to 100 yrs). We see that the quadrupole 
approximation becomes less accurate when we increase the 
semi-major axis ratio a. This could be due to truncation 
of the Hamiltonian at order q^, or it coul d be due to in- 
accur acies of the first order secular theory jGiuppone et al.l 
[ioil). These results do not vary much with other choices 
for the masses mo and mi, as long as mt = mo + mi is kept 
constant. 



3 ORBITAL PRECESSION IN RADIAL 
VELOCITY DATA 

3.1 Measuring precession rates 

The radial velocity of a star, m2, with a close binary system, 
mi, — mo + mi , is approximately 



Vr = VrO + Vri 



(27) 



where Vrat are short per iod perturbation terms obtained in 
iMorais fc Correial (|201ll '). 

(28) 
(29) 



VrO = K{cOs{f2 + 0^2) + 62 COs(tJ2)) 

with 

71-2 02 mt 



K = 



- e| "^2 -I- mb 



sin/2 , 



and U12 changes linearly with time (Sect. 2.3), i.e. 

U2 — lJJ20 +0J2t . (30) 

The precession of the star's orbit can be inferred from 
the radial velocity curve mostly due to the term oc 62 in 
Eq. (|28|) . Typically, the observation timespan, tob3, is much 
shorter than the precession cycle, thus if [J20 7^ 0, 180° 




time (years) 

Figure 3. Evolution of angle uJ2, showing secular drift and short- 
period oscillations, obtained for simulations as in Sect. 4.1 but 
without observational errors; case I (top panel) and case II (bot- 
tom panel). 



K 62 cos(cj2) ^ K €2 [cos(tt;2o) — sin(cJ2o) dj2 t] (31) 

with t < tabs- 

Hence, Eq. (|28|l is approximately a Keplerian radial ve- 
locity curve whose amplitude has a linear drift which is at 
most {0J20 = 90°, 270°) 



K 62 Lb2 tob 



(32) 



If UJ20 ~ 0, 180°, the radial velocity curve's amplitude has a 
quadratic drift 7^62 (tj2 tobs)^/2. 

In order to measure ijj2 with accuracy, two conditions 
must be met. First, the drift (Eq. (O) must be larger than 
the observation's precision. Second, the observation times- 
pan, tabs, must be a few outer binary's periods so that we 
can distinguish the sec ular drift, ijj2toba , from short period 
oscillations, Alj (Moral s fc Correiall201lh . i.e. we must have 
ijj2 tabs 3> Ao; (see Fig. 3) . These two conditions help us 
predict when can we measure accurately the outer binary's 
precession rate. However, in practice we estimate d}2 by fit- 
ting a precessing Keplerian orbit (Eqs. [28] and I30p to the 
radial velocity data. 



3.2 Estimating inner binary parameters 

We saw that (Eqs. ^ and ((ggj) 

u;2 oc x(l — x)aiA . 



(33) 
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where x = mi/nit, mt — mo + mi and A is a function of 6 
(cf. Sect. 2.3). 

Therefore, we can use the measured precession rate, ijj2, 
to estimate the parameters of an hidden binary (i.e. mi and 
ai) through the quantity x{l — x)af. However, since A is 
a function oi 6 = cosi which is unknown, we have to make 
some assumptions. From Fig. 2 in Sect. 2.3 we see that pro- 
grade precession is faster when i — and that, if the inner 
binary is at the Kozai stationary solution, retrograde preces- 
sion occurs if i > 45° and it is faster when i — 90° . As orbits 
with high relative inclination, i, will reach values of ei near 
unity, these are likely to become unstable or undergo tidal 
evolution, hence we set a maximum value of i = 60° . There- 
fore, if tj2 > we assume that i — and ei = 0; if a;2 < 
we assume i — 60° and ei — 0.76. These assumptions imply 
maximum precession ratefl (in absolute value), thus they 
give us minimum estimates for the parameter x{l — x) a\. 

In order to estimate the hidden inner binary compo- 
nent's mass, mi, we must provi de an estimate for t he in - 
ner binary's semi-major axis, a\ . iHolman k, WiegertI (| 19991 ) 
measured the size of stability regions around binary star 
components. They assume a massless particle orbiting in the 
binary system's plane and find that this is stable if a\ < ac 
with 

— = (0.464 ± 0.006) + (-0.380 ± 0.010) 

0,2 

-f(-0.631 ± 0.034) £2 + (0.586 ± 0.061) ^ 62 

-f- (0.150 ± 0.041) el + (-0.198 ± 0.074) /lei (34) 

where a2 and 62 are, respectively, outer binary's semi-major 
axis and eccentricity, and mass parameter — m2 / (m2 + 
mb). 



4 EXAMPLES 

To test our model, we performed numerical integrations of 
the equation of motion of hierarchical triple systems com- 
posed of a star, m2, and a close by binary, mt- We chose 
masses m2 — Mq and mb = 0.5 Mq, and semi- major 
axes a2 = 3 AU and ai = 0.3 AU, which imply periods 
T2 = 4.24 yr for the outer binary and Ti = 84.9 day for the 
inner binary, and semi- major axis ratio a — ai/a2 — 0.1. 
The initial angles were h = 90°, fli ^ ^2 = 0, lji = 90°, 
u}2 = 20°. We computed the radial velocity of the star m2, 
and simulated observational data points for a timespan, tabs, 
and a certain precision limit. We then applied the traditional 
techniques used in radial velocity data analysis. 



4.1 Measuring precession rates 

To test in which circumstances we are able to measure the 
outer binary's precession rate we set the outer binary on an 
eccentric orbit (62 = 0.2). 

Case I is a coplanar triple system (i = i.e. Ii = I2 = 
90°) where the inner binary has masses mo = 0.42 A/0 and 
mi = 0.08 Mq, and a nearly circular orbit (ei = 0.01). In 



* However, if ei 7^ 0, the precession rate at i = increases by a 
factor 1 -(- 3 /2 with respect to ei = 0. 



Case I 


fit 


A 


B 


prcc (m/s) 




5.425 


1.085 




(1) 


1544.12 ± 0.05 


1 544 15 + 001 




(0) 


1544.15 ±0.05 


1544.18 ± 0.01 


K (m/s) 


(1) 


7172.9 ± 0.7 


7173.86 ± 0.13 




(0) 


7172.6 ±0.7 


7173.62 ± 0.13 


e 


(1) 


0.19879 ±0.00009 


0.19873 ±0.00002 




(0) 


0.19875 ±0.00009 


0.19870 ±0.00002 


uj (Vyr) 


(1) 


0.115 ±0.016 


0.0937 ±0.0032 




(1) 


1.569 


1.660 




(0) 


1.690 


3.097 


rms (m/s) 


(1) 


9.1850 


1.9154 




(0) 


9.8849 


3.4024 



Table 1. Fits to Case I (mo = 0.42 M©, mi = 0.08 M©, i = 0, 
ei = 0.01) with toba = 8 yr and different precisions. 



Table 1 we show the results of fitting a fixed Keplerian or- 
bit (fit 0) or a precessing Keplerian orbit (fit 1) to the data 
time series with 128 points over tabs ~ 8 yr, at precisions of 
about 5 m/s (A) and 1 m/s (B), respectively. From Fig. 3 
(top) we see that since tabs = 8 yr, Auj > ij2 toba- The maxi- 
mum radial velocity drift over tabs is 18.6 m/s. At precision 
5.425 m/s (A) the observation error is 29% of the maximum 
drift. Therefore, fit (1) is only slightly better than fit (0). 
However, at precision 1.085 m/s (B) the observation error is 
only 6% of the maximum drift, thus fit (1) is clearly better 
than fit (0). The theoretical value (quadrupole approxima- 
tion) for the precession rate is uj2 = 0.093° /yr while the true 
value (simulations up to 100 yrs) is UJ2 ~ 0.089° /yr. 

Case II has an inner binary with masses mo = 0.35 Mq 
and mi = 0.15 Mq, with i = 60° (i.e. h = 30° and h = 90°) 
and ei — 0.76 (Kozai stationary solution). In Table 2 we 
show the results of fitting a fixed Keplerian orbit (fit 0) 
or a precessing Keplerian orbit (fit 1) to data time series 
with 99 points over tabs ~ 6 yr (C), and to data time series 
with 154 points over tabs ~ 12 yr (D), both at precision 
of about 5 m/s. The maximum radial velocity drifts over 
6 yr and 12 yr are, respectively, 31 m/s and 62 m/s, which 
correspond to observation errors of, respectively, 16% and 
8% of the maximum drift. From Fig. 3 (bottom) we see 
that when toba = 6 yr, Aui « 6j2 tabs while when toba = 
12 yr, Acj <C ijJ2toba- Therefore, when toba = 6 yr (C), fit 
(1) is slightly better than fit (0) but when toba ~ 12 yr 
(D), fit (1) is clearly better than fit (0). The theoretical 
value (quadrupole approximation) for the precession rate is 
ijj2 = — 0.272° /yr while the true value (simulations up to 
100 yrs) is [1)2 = -0.238° /yr. 

4.2 Estimating inner binary parameters 

We simulated triple systems as described above with 
the observed star on an eccentric orbit with 62 = 0.2 
(outer binary) around an inner binary with masses mo — 
0.42 Mq and mi = 0.08 Mq, and semi-major axis ai = 
0.3 AU. The values of the relative inclination were i = 
0°, 20°, 40°, 50°, 55°, 60°. When i < 40° the inner binary 
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Case II 


fit 


c 


D 


tabs (yr) 




6 


12 






1 567 74 + 07 


1 567 75 + 03 




(0) 


1567.44 ± 0.07 


1567.70 ± 0.03 


K (m/s) 




7135.81 ± 0.89 


7135.61 ± 0.60 




(0) 


7139.36 ±0.84 


7134.50 ± 0.60 


e 


(1) 


0.20495 ±0.00011 


0.20504 ± 0.00009 




(0) 


0.20491 ±0.00011 


0.20558 ± 0.00009 


^ (7y) 


(1) 


-0.256 ±0.021 


-0.2606 ± 0.0086 




(1) 


1.656 


1.628 




(0) 


2.077 


2.984 


Tvas (m/s) 


(1) 


9.4952 


9.5135 




(0) 


11.6178 


16.1299 



Table 2. Fits to Case II (mo = 0.35 M©, mi = 0.15 Mq, i = 60°, 
ei = 0.76) with precision 5.4 m/s and different tofts- 
had a nearly circular orbit (ei = 0.01) and when i > 40° it 
had an eccentric orbit near the Kozai stationary solutior0. 
In all cases tots ~ 8 yr, and the precision limit was about 
1 m/s. 

In Table 3 we present the precession rates (theoreti- 
cal and measured in the simulations) and the ratio between 



of fit (1) and fit (0) which measures the goodness of 
fit (1) with respect to fit (0). As explained in Sect. 3.2, we 
obtain minimum estimates for the inner binary parameter, 
x{l — x)al with X = mi/mj, (Eg. I33|) . assuming i — and 
ei = when tj2 > 0, or i = 60° and ei = 0.76 when (jj2 < 
(cf. Fig. 2). We can then obtain minimum estimates for 
mi assuming ai = Oc = 0.49 AU which is the maximum 
size of stable orbits (massless particle mi <^ mt) around 
mo « mb in the coplanar case (see Eq. (I34p with 62 — 0.2 
and /i = 0.67). These estimates are all realistic (minimum 
mass of the hidden inner binary companion between 6 Mj 
and 27 Mj) hence can be used as input parameters for a 
N-body fit which can provide best-choice values for mi and 
ai. 



4.3 Can precession mimic a planet? 

Here, we repe a t the qu estion already made in 
iMorais fc Correial (|2008l . I2OIII ). If we do not know about 
the inner binary's presence because one of its components is 
unresol ved, can the binar y 's effect be mistaken as a planet? 

In ' Morals fc Correial (|2008l . 120111 ) we showed that the 
residuals leftover from fitting a fixed Keplerian orbit to the 
outer binary (observed star's orbit around the inner binary's 
centre of mass) contained additional periodic signals that 
could be mistaken by planets. Here, we show an example 
of similar behavior obtained from the previous simulation 
(case II with tabs = 6 yr) in Fig. 4 (top). The periodograrrff] 



i (°) 


ei 


dJ2 (t) 


1^2 (s) 


ratio 


mi (Mq) 





0.01 


0.0928 


+0.094 
±0.003 


0.54 


0.027 


20 


0.01 


0.0765 


+0.078 
±0.003 


0.60 


0.022 


40 


0.01 


0.0353 


+0.040 
±0.003 


0.85 


0.011 


50 


0.56 


-0.0473 


-0.040 
±0.003 


0.89 


0.006 


55 


0.67 


-0.1065 


-0.099 
±0.003 


0.67 


0.015 


60 


0.76 


-0.1740 


-0.167 
±0.003 


0.50 


0.025 



^/l — (5/3) cos^ i and 



^ The Kozai stationary solution has e\ 
cji = ±90°. 

^ We c ompute Generalized Lomb-Sc argle periodograms as de- 
fined in IZechmeister fc KursteJ | |2009|) 



Table 3. Hierarchical triples system (a = 0.1) with inner binary 
on nearly circular orbit (if i < 40°) or at Kozai stationary solution 
(if i > 40°): comparison between outer binary's theoretical (t) and 
observed (s) precession rates; y/x^ ratio between fit (1) and fit 
(0) to outer binary's orbit; minimum mass of hidden inner binary 
component, mi. 



has an obvious peak at 606 day which is nearly commensu- 
rate (ratio 2/5) with the outer binary's period. However, this 
peak disappeared when we fitted a precessing Keplerian or- 
bit to the outer binary (Fig. 4: middle), and it was no longer 
prominent when we increased the observation timespan to 
tobs = 12 yr (Fig. 4: bottom). 

In|M orais fc Correial l|201lf ) we studied the short term 
effect of a binary system on a nearby star in the case of ec- 
centric and inclined orbits. We saw that, when the observed 
star's orbit (outer binary) was eccentric, the radial velocity 
was composed of a main Keplerian term that described the 
star's motion around the inner binary's centre of mass, and 
short-period term s (obtained by integratin g with respect to 
time Eq. (34) in iMorais fc Correial HoiJ)). In particular, 
some short-period terms appeared at harmonics of th e outer 
binary's frequency, 712. However, in lMorais fc Correial (|201ll ) 
we considered that tabs was only a fraction of the outer bi- 
nary's period. In this situation, the outer binary's orbit was 
not well constrained and these harmonics were incorporated 
into the main Keplerian term. Now, as tabs covers a few outer 
binary's orbits, these harmonics will appear in the residu- 
als leftover from fitting a fixed Keplerian curve to the outer 
binary. However, the observed star's orbit (outer binary) is, 
in fact, precessing and the precession rate, aJ2, is approxi- 
mately constant in the long-term but has short-term oscil- 
lations (see Fig. 3). If tabs covers only a few outer binary's 
orbits, these short term oscillations will cause mixing up of 
the frequencies thus the signals do not exactly coincide with 
harmonics of n2. In particular, combinations of these har- 
monics can appear (Fig. 4: top). When we fit a precessing 
Keplerian orbit to the outer binary, as tabs is short compared 
to the precessional period, the signals at or nearby harmon- 
ics of 712 are incorporated into the precessing Keplerian orbit 
(Fig. 4: middle). As tabs increases, the short period oscilla- 
tions become negligible with respect to the secular terms 
(Fig. 4: bottom). 

If the observed star has a circular orbit then there is 
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Figure 4. Periodogram of residuals after fit (0) and fit (1) to 
radial velocity data from simulation in Table 2 with t^i,^ = 6 yr 
(top and middle, respectively) and after fit (0) to simulation in 
Table 2 with ti^i,^ = 12 yr (bottom). 



no pericentre precession. When the inner binary's orbit is 
also circular but inclined {i < 40°) with respect to the outer 
binarjQ, the star's radial velocity is composed of a main Ke- 
plerian term (circular orbit with frequency ^2) and short 
period term s (obtained by int e gratin g with respect to time 
Eq. (21) in iMorais fc Correial (|201ll )). In particular, there 
are signals at frequencies n2 and 3n2. The term with fre- 
quency 722 is simply incorporated into the main Keplerian 
fit. However, the term with frequency 3n2 can be mistaken 



^ Due to the Kozai effect, when i > 40° the inner binary's orbit 
cannot remain circular. 




Figure 5. Periodogram of residuals after fit (0) to radial velocity 
data from simulation of circular non-coplanar triple star system 
composed of an observed star and a companion unresolved binary. 
Peaks at 46 day and 505 day can be mistaken by planets. 



as a planet at the 3/1 mean motion resonance with a com- 
panion "star" of mass irib — mo + mi. 

We simulated a triple system, as explained at the begin- 
ning of Sect. 4, with mi — m2 = 0.25 Mq , 62 = 0, ei = and 
i — 30°. We generated radial velocity data with 154 points 
over tabs = 12 yr, at precision 0.543 m/s. In Fig. 5 we show 
a periodogram of the residuals leftover after fitting a Keple- 
rian orbit to the outer binary. As expected, we see peaks at 
505 day (frequency 3 712, harmonic of n2 ) and 46 day (fre- 
quenc y 2 Til — 3 712 , short-period term as in lMorais fc Correial 
l|2008l )) with amplitudes 1.7 m/s and 1.2 m/s, respectively. 
These can be mistaken as planets. 



5 A PLANET IN I/-OCTANTIS? 

The system i/-Octantis is a close single-line spectroscopic bi- 
nary. |Kamm et all (|2009l ) published radial velocity data con- 
sisting of 221 points covering a timespan tots = 1862 day, 
with precision around 5 m/s (inferred from the published 
observation errors). Combi ning the radia l veloc ity data with 
astrometric measurements, iRamm et al.l l|2009l ) derived im- 
proved parameters for the i/-Octantis binary (Table 4). The 
residuals leftover from fitting a Ke plerian orbit to ty-O ctantis 
show an additional signal, which iRamm et al.l l|2009l) iden- 
tify as a planet with minimum mass, 2.5 Mj, and semi-major 
axis, 1.2 AU, orbiting the primary star (cf. Table 4). More- 
over, the planet's period, 417 day, is nearly commensurate 
(ratio 2/5) with i/-Octantis binary system's period. 

A planet about half way between the primary and sec- 
ondary stars (semi- major axis ratio a = .47) i s unexpected. 
In fact, according to lHolman fc WiegertI (Il999l ) the stability 
limit for coplanar prograde orbits around the primary star 
is only = 0.6 AU {Eq.^. Nevertheless. lEberle fc Cuntj 
(|2010l ) propose that such planet (op = 1.2 AU) can be stable 
for at least 10 million years on a retrograde coplanar orbit. 

In Table 5 we pr esent the results of fitting the radial 
velocity data from iRamm et al] (|2009) with a fixed Keple- 
rian orbit (fit 0) and a precessing Keplerian orbit (fit 1). 
We see that fit (1) with LO2 = -0.86 ± 0.02°/yr- is better 
than fit (0), although the difference is not yet very signif- 
icant. However, as seen previously in Sect. 4.1, this could 
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v-Octantis A 



rriA = (1.4±0.3)Mo 



i^-Octantis B + planet 



T2 = (1050.11 ±0.13) day 


Tp = (417 ± 4) day 


K2 = (7032.3 ± 2.6) m/s 


Kp = (51.8 ± 1.6) m/s 


62 = 0.2359 ±0.0003 


ep = 0.123 ±0.037 


UJ2 = (75.05 ± 0.075)° 


uii = (260 ± 21)° 


h = (70.8 ±0.9)° 


Ip=? 


Ha = (87 ±1)° 


np =? 


ms = (O.5±O.1)M0 


TUp sin/p = 2.5 M J 


a2 = (2.55 ±0.13) AU 


ap = (1.2 ±0.1) AU 




4.2 


rms = 19 


m/s 



Table 4. Fitted p arameters for 
llRamm et al.ll2009l) . 



'-Octantis and possible planet 



be due to the short observation timespan {tabs = 1.77 Ta). 
In Fig. 6 we show the periodogram of the residuals after fit 

(0) and fit (1). After fit (0) there is a prominent signal at 
417 day, while after fit (1) this signal is still present but it is 
no longer dominant and seems to be within the noise level. 
This is similar to the behavior described in Sect. 4.3. From 
Table 5 we see that although fit (1) is better than fit (0), the 
fit with a planet at 417 day is currently better. As explained 
above, this could be due to the short-observation timespan, 
or even due to the particular sampling of the radial velocity 
data. Moreover, the fit with the planet introduces additional 
five free parameters while fit (1) introduces only one more 
free parameter (0)2) which could also help explain why the 
fit with the planet seems better than fit (1). 

We saw that retrograde precession occurs if the sec- 
ondary star in i^-Octantis is in turn a binary system in- 
clined more than 45° with respect to the main binary's orbit 
(Fig. 2). In this scenario, we saw that a periodogram of the 
residuals leftover from fitting a Keplerian orbit to the main 
binary could exhibit peaks that might be mistaken as plan- 
ets. We saw that these peaks appeared close to harmonics of 
the main binary's orbital frequency. In particular, we showed 
an example (Sect. 4.3 and Fig. 4) where the fake planet's 
period is nearly commensurate (ratio 2/5) with the main bi- 
nary's period which is exactly what happens in !/-Octantis. 
Therefore, we propose that a hidden binary system could 
mimic a planet sim ilar to the one reported in i/-Octantis by 
iRamm et al] (|2009l ). 

In order to estimate inner binary parameters that lead 
to retrograde precession of the outer binary at a rate 0^2 ~ 
-0.86° /yr, we set ai = = 0.35 AU (Eq. [Ml), « = 60° 
and ei = 0.76 (i.e. the inner binary is at the Kozai station- 
ary solution). Replacing these in Eqs. (25) and (26) with 
02 = 2.55 AU, 7712 = 1.4 Mq and mt, = 0.5 Mq, we obtain 
estimates for the inner binary's masses of mi = 0.23 Mq and 
mo = 0.27 M0. In Table 5 we show the results of fitting a 
precessing Keplerian orbit (fit 1) and a fixed Keplerian orbit 
(fit 0) to such triple system (simulation I). We see that fit 

(1) is better than fit (0) and is comparable with the fit of a 
planet at 495 day. 

We can also estimate inner binary parameters by per- 




Figure 6. Fourier analysis of residuals after fit (0) and fit (1) to 
i^-Octantis radial velocity data (top and bottom, respectively). 



forming N-body fits to i^-Octantis radial velocity data. We 
assumed m2 = 1.4 Mq, and fixed the outer binary's orbit 
at I2 = 70.8° and ^^2 ~ 87° which are the parameters 
inferred from the spectroscopic-astrometric solution (Table 
4). We present the best fit solution in Table 6. Retrograde 
precession occurs because i = 122°. Moreover, the best bi- 
nary N-body fit (Table 6) is comparable to the Keplerian 
planet fit (Table 4). Although this solution is unstable be- 
cause ai ~ 1.5 Oc, we can obtain "equivalent" stable config- 
urations by reducing ai, while maintaining nib ~ 0.538 
constant, and increasing the ratio mi/mb, so that the Ke- 
plerian term and the quadrupole interaction term are both 
kept constant. 

We performed a simulation of such a stable configu- 
ration with tti = flc = 0.35 AU, mi = 0.109 Af© and 
mo = 0.429 Mq . In Table 5 we show the results of fitting 
a precessing Keplerian orbit (fit 1) and a fixed Keplerian 
orbit (fit 0) to such triple system (simulation II). We see 
that fit (1) is better than fit (0) but fit (1) is better than 
the fit of a planet at 452 day. The values of the precession 



rate after fit (1), and both ^^'^ residuals {rms) after 

fitting the planet's orbit are almost equal to those obtained 
for the real system j/-Octantis (cf. Table 4). 

From Table 5, as described above, we see that for sim- 
ulated data the fit of a precessing Keplerian orbit is com- 
parable (simulation I) or better (simulation II) than the fit 
of a planet, while in the real case (iz-Octantis) the fit of a 
planet is currently better than the fit of a precessing Kep- 
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fit 


//-Octantis 


simulation I 


simulation II 




(1) 


1050.46 ± 0.03 


1114.63 ± 0.04 


1093.44 ± 0.03 


T (day) 


(0) 


1050.11 ± 0.03 


1114.32 ±0.03 


1093.04 ± 0.03 




(0)+pl 


417 ±1 


495 ±2 


452 ±2 




(1) 


7044.24 ± 0.60 


6895.7 ±0.5 


6970.6 ±0.6 


K (m/s) 


(0) 


7032.27 ± 0.68 


6889.18 ±0.69 


6961.7 ±0.6 




(0)+pl 


51.83 ± 0.53 


18.65 ± 0.69 


lb. 2. ± 0.7 




(1) 


0.23553 ± 0.00007 


0.24946 ±0.00008 


0.25388 ± 0.00008 


e 


(0) 


0.23589 ± 0.00009 


0.24767 ±0.00011 


0.2524 ±0.0001 






124 + 010 


51 + 03 


33 + 02 


(7yr) 


(1) 


-0.860 ±0.017 


-0.500 ±0.017 


-0.810 ±0.016 




(1) 


7.3 


2.374 


3.6 




(0) 


8.1 


3.098 


5.0 




(0)+pl 


4.4 


2.364 


4.3 




(1) 


36.3 


13.54 


19.3 


rms (m/s) 


(0) 


39.1 


16.92 


26.7 




(0)+pl 


22.8 


13.30 


22.8 



Table 5. Fits to !^-Octantis (left) and simulations I and II (right). The radial velocity data in I 
and II was generated at the 221 observational data points for i/-Octantis, covering toi,^ = 5.1 y 
with precision about 5 m/s. 



!^-Octantis A 



mA = 1.4 Mq 



f-Octantis B = hidden binary system 



T2 = (1078 ± 1) day Ti = (189.1 ± 1.4) day 

K2 = (7010 ± 3) m/s Ki = (2812 ± 76) m/s 
62 = 0.2504 ± 0.0003 ei = 0.67 ± 0.03 

UJ2 = (72.63 ± 0.13)° uii = (23.25 ± 1.38)° 



h = 


(70.8 ±0.9)° 


h = (63.4 ±2.7)° 




= (87 ± 1)° 


f^l = (232 ± 1)° 


mo 


= 0.496 A/0 


mi = 0.042 Mq 


112 


= 2.565 AU 


ai = 0.524 AU 






= 4.9 




rms 


= 25.9 m/s 



Table 6. Best N-body fit for i/-Octantis assuming hierarchical 
triple star system. 



lerian orbit. However, we stress that there are many more 
combinations of parameters (mo, mi, i and ei) that can 
cause a precession rate uj2 = — 0.86°/yr, assuming that it is 
well constrained. Moreover, we expect our model to explain 
better the synthetic data generated with 3-body simulations 
than real data (!/-Octantis) where we could have other plan- 
ets or even stellar variability. 

A planet around the primary star in //-Octantis can also 
cause precession of the main binary's orbit. However, our 
simulation s show that coplanar re trograde planet orbits, as 
reported in lEberle fc Ciint3 l|20ld ). cause slow prograde pre- 
cession of the main binary's orbit at a rate 0.04° /yr. This 
is also what we expect from our quadrupole order theory 
(Eqs. 1251 126() although we do not expect it to be accurate 
at semi-major axis ratio a = 0.47. We saw (Sect. 2.3) that 



in order to have retrograde precession we would need the 
planet's orbit to be inclined more than 45° with respect 
to the z/-Octantis binary. In our numerical integrations we 
could not find (although we did not do an exhaustive search) 
stable planet orbits at semi-major axis ratio a — 0.47 and 
with such high inclination with respect to the main binary. 



6 CONCLUSION 

We studied the effect of a binary system on a nearby 
star's motion. This is a comp lement of our previous work 
iMorais fc Correial l2008l . l201ll ) where we assumed that we 
had observations for a fraction of the star's orbit around the 
binary's centre of mass. Here, we assumed that we had ob- 
servations for a few orbits of the star around the binary's 
centre of mass. We saw that, in this case, the secular effect 
of the binary dominates over the short-term effects. 

We developed a secular theory which was based on a 
quadrupole expansion of the Hamiltonian. This is accurate 
for hierarchical triple systems composed of an inner binary, 
and a star that moves around this inner binary's centre of 
mass on a wider orbit which we called the outer binary. 

We derived an expression for the outer binary's pre- 
cession rate and showed that it is approximately constant. 
Therefore, the star's radial velocity can be modeled as a 
modified Keplerian radial velocity curve with slowly drift- 
ing amplitude. We then showed how we can measure the 
outer binary's precession rate by fitting a precessing Keple- 
rian orbit to the radial velocity data. We also showed how 
we can estimate inner binary parameters from the measured 
precession rate. 

We saw that, if we are unaware of the inner binary's 
existence and simply fit a non-precessing Keplerian orbit to 
the radial velocity data, a periodogram of the residuals will 
show peaks at or nearby harmonics of the outer binary's 
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period which can be mistaken as planets. However, if we fit 
a precessing Keplerian orbit to the radial velocity data, these 
signals are no longer prominent in the leftover residuals. We 
conclude that detecting precession in the radial velocity data 
of a star within a binary system may be an indication that 
there is an unresolved third star. 

We discussed the case of z/-Octantis which is a close 
binary system (2.55 AU) composed of a K-type star (i/- 
Octantis A) and a fainter companion (//-Octantis B). Radial 
velocity data analysis showed a signal at 417 day which was 
identified as a planet at 1.2 AU of iz-Octantis A. However, we 
showed that the radial velocity data currently implied retro- 
grade precession of about — 0.86°/yr for this binary system. 
We suggested that this may indicate that z/-Octantis B is 
actually a double star which could explain a signal similar 
to that previously associated with a planet. At the moment 
we cannot yet decide that the reported planet of (/-Octantis 
A is simply an artifact caused by i/-Octantis B being a dou- 
ble star. In order to distinguish between the two hypoth- 
esis (planet or double star), more radial velocity data for 
i/-Octantis is needed, so that we can better constrain the 
main binary's precession rate. Moreover, the planet hypoth- 
esis could be compatible with retrograde precession of the 
i^-Octantis binary if we could prove the existence of sta- 
ble orbits around i^-Octantis A, with semi-major axis ratio 
a — 0.47 and inclined more than 45° with respect to the 
I'-Octantis binarjjfj. 
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* As we were in the process of submitting the final version of 
this article, we became aware of a stability study done by K. 
Gozdziewski which confirms that stable planetary solutions com- 
patible with i^-Octantis radial velocity data are unlikely. 



